The study of random vorticity in quantum fluids through interference fluctuations 
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We study the vortex dynamics of a quantum degenerate Bose gas through the intensity fluctuations 
of the interference from particles extracted at two different positions. It is shown numerically 
with classical field simulations that an interacting Bose gas with proliferating vortices exhibits long 
correlation times for these intensity fluctuations. This behavior is contrasted with the case of a 
noninteracting gas, that we describe analytically, and with the case of a well condensed Bose gas 
without vortices. We discuss the observability of our predictions in quantum fluids of exciton- 
polaritons. 
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I. INTRODUCTION 

Quantized vortices play a central role in the thermal 
and dynamical properties of quantum degenerate Bose 
gases. Vortex lattices in rotating Bose-Einstein con- 
densates [1 , spontaneous vortices above the Berezinskii- 
Kosterlitz-Thouless (BKT) transition [2 , discrete vortic- 
ity in turbulent superfluids [3 and Kibble-Zurek (KZ) 
vortices [1] are prominent examples of the importance of 
quantized vortices. 

Their experimental visualization relies either on the 
associated density minimum or the singularity in the 
phase profile. Where the density is usually the more 
easily accessible observable, only a measurement of the 
phase gives unambiguous proof of the quantized circu- 
lation around the vortex core. Vortices in both Bose- 
Einstein condensates of ultracold atomic gases and of 
exciton polaritons in semiconductor microcavities have 
been detected in this way. To this end, the particles 
coming from different regions in the same condensate (of 
exciton-polaritons [5 ) or from independent condensates 
(of atomic gases [2, 6 ) were interfered. A dislocation in 
the interference pattern proves the presence of a quan- 
tized vortex. With this technique, vortices in rotating 
atomic condensates have been evidenced [6^, as well as 
spontaneous vortices above the transition temperature 
of a 2D atomic Bose gas [2]. In exciton-polariton con- 
densates, quantized vortices induced by the interplay of 
pumping and decay with a disordered potential landscape 
were discovered 0. 

Regarding their observation, the distinction between 
deterministic and spontaneous vortices is essential. The 
former occur on positions that are determined by the sys- 
tem geometry - examples are the vortex lattices of ro- 
tating condensates or the vortices of exciton-polaritons 
pinned to disorder. Their predictable positions makes 
that there are no restrictions on the time it takes to 
record the interferrogram. For the observation of the 
spontaneous quantized vortices on the other hand, it is 
essential that the experiment can be performed on a time 
scale that is fast with respect to the vortex motion. For 
experiments with ultracold atoms, where the intrinsic 



time scales are of the order of ms, this poses no prob- 
lem. The spontaneous BKT [2 and KZ [4 vortices have 
indeed been successfully detected. For exciton-polariton 
condensates [7|, that are a billion times faster (typical 
time scale is 1 ps), the required detection speed poses a 
technological challenge that is yet to be overcome. The 
time it takes to record an interferrogram is much longer 
than the time scale on which the vortices move. As a con- 
sequence, after the temporal averaging, no phase defects 
can be observed, but only a decreased spatial coherence. 
This difficulty has made spontaneous hydrodynamic vor- 
tices under cw excitation [8^ and possible KZ vortices [9 
elusive and forms a limitation for the experimental study 
of quantum turbulence, suggested in Ref. [10]. 

While no image of the polariton interferrograms can 
be recorded on a ps time scale, it is possible to study its 
intensity fluctuations with a standard Hanburry Brown 
and Twiss setup [TP. We will propose in this paper to 
study these intensity fluctuations in order to characterize 
the vorticity in the fluid. The idea is presented in Sec. 
im Where the intensity fluctuations can be computed an- 
alytically in the noninteracting case, numerical classical 
field simulations are presented in Sec. [TIT] to investigate 
the interacting case. The effect of a finite particle life 
time, relevant for experiments with exciton-polaritons is 
addressed as well. We draw conclusions in Sec. IIVI 



II. INTERFERENCE INTENSITY 
FLUCTUATIONS 

For a schematic understanding of our proposal, we have 
sketched in Fig. [l] an interferrogram of a condensate con- 
taining a vortex with a displaced copy of itself. Due to a 
TT phase shift between the two interfering beams, a non 
vanishing intensity is visible only between the positions 
of the vortices, that are marked with a full circle (vortex 
position on the condensate image) and an empty circle 
(vortex position on the displaced image) . Panel (c) shows 
a separation of the interfered images a that is smaller 
than the vortex core size ^. The maximal density in the 
interferrogram is much smaller in this case than in panel 
(d), where the separation a is larger than the vortex size. 
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FIG. 1: Sketch of the density that results from interfering a 
condensate containing a vortex (a) with a displaced copy (b) , 
with a phase shift of tt between the two. The interferrogram 
shows a density maximum between the vortex cores, whose 
size is determined by the displacement a. Panels (c) and (d) 
show two distances a, that are smaller and larger than the 
vortex core size respectively. 



The reason is simply that for small separation, construc- 
tive interference occurs only in the vicinity of the vortex 
where the density is suppressed. 

It is expected that vortices that move ballistically 
should give rise to correlations in the density fluctua- 
tions over a time that scales as r^ ~ a. In the case of 
diffusive vortex motion on the other hand, one rather ex- 
pects td ~ cl'^' Below, we will show analytically that 
this scaling constitutes a qualitative difference from the 
case of the noninteracting gas where the correlation time 
depends much weaker on a. The scaling Tc ~ a^a^ in 
the presence of vortices will be addressed with numerical 
simulations in a classical field model. 

When the light extracted at the positions x and x + a 
is interfered, the intensity incident on the detector is 
^61 = V^(x) + e^^?/;(x + a). When 6 is chosen so to min- 
imize the amplitude of ^, the intensity on the detector 
is the smallest and the relative fluctuation should be the 
largest. We restrict to this case in the following of the 
paper. For ordinary condensation in a k = state, this 
corresponds to 6> = tt, but in the general case, a phase 
that suits the condensate wave function should be chosen. 
In typical experiments with polariton condensates [T2|, 
this means that the density fluctuations are measured at 
a minimum of the interference fringes. The second order 
correlation function 

5f (a, t, r) ^ {^Umiit + T)*,(t + r)vl/.(i)) (1) 

quantifies the correlations between a density fluctuation 
at time t and one at time t + r. 

When the fluctuations on top of the condensate are 
Gaussian, the density fluctuations of the interferrogram 
are Gaussian as well. We then find in the steady state a 



bunching of the fluctuations: g\ ^(a, t = 0) = 2g\ ^(a, t = 
oo). When the different momentum components are as- 
sumed to have no phase relation, the decay time of the 
interference intensity fluctuation can be computed ana- 
lytically by using Wick's theorem (see appendix) 
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where e(/c) is the single particle dispersion. The nor- 
malization constant A/" is determined such that |/(r = 
0)1 = 1. Eq. ^ expresses a simple relation between the 
density fluctuations in the interferrogram and the mo- 
mentum distribution of the fluid. 

The validity of Eq. Q depends on the absence of cor- 
relations between the different momentum components. 
This condition is satisfied for the non-interacting Bose 
gas. When interactions are turned on, the momentum 
components of the gas become correlated, affecting the 
fluctuations of the interferrogram. We will study this 
case numerically in the next Section. 



III. CLASSICAL FIELD SIMULATIONS 

We will investigate the weakly interacting Bose gas nu- 
merically based on a classical field model[13 . In this 
method [14 , the quantum field is replaced by a classical 
field whose dynamics is described by the Gross-Pitaevskii 
equation (GPE) 
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Here, we have used a quadratic disperion with mass m 
and g is the strength of the contact interaction. The 
evolution of the classical field with the GPE gives a rea- 
sonably accurate description of the thermal bose gas, in- 
cluding the BKT transition [15 . In this work, we do not 
to focus on the spontaneous vortices of the BKT type. 
Actually in the the weakly interacting Bose gas, no well 
defined spontaneous vortices exist |T6], due to the large 
magnitude of density fluctuations close to the transition 
temperature. 

It is still possible to construct a weakly interacting 
quasi-condensate with well defined vortices by a suit- 
able initial condition. Experimentally such a state can 
be the result of a sudden temperature quench of a gas 
of ultracold atoms [T7], that produces vortices of the 
Kibble-Zurek type [4 . In the polariton case, the ini- 
tial state can be chosen at will, because under reso- 
nant excitation, it is directly imprinted by the excita- 
tion laser. A different approach is a quantum turbu- 
lent state where vorticity forms when the bose gas is 
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FIG. 2: Phase (a) and density (a) profile of a snap shot of the classical field that is evolved with the GPE for t — 50ps from 
a Gaussian momentum distribution (width 0.5/im~^) without interactions (I), with interaction strength gn — 0.5meV (II) and 
from a condensed state (96 percent condensate) with Gaussian occupation of higher momentum states (III). Panels (c) show 
the corresponding Qj (normalized to its value at zero delay) as a function of the delay r for several values of the separation a. 
The panels (d) show the same, but for a rescaled time variable r/a. The inset in panel lid) shows ^j as a function of r /a^ . 
Other parameters: effective mass h/m — l/xm^meV, interaction energy gn — O.^meV . Simulations were performed on a system 
of size 100x100 /xm and 128x128 grid points. 
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driven externally [TDl [TH]. In order to study a quasi- 
condensate with vortices, we choose a simple Gaussian 
initial condition in momentum space with random phases 
for each state '^(/c) = exp(— /c^/Zcq + irjik))^ where k^ de- 
termines the width of the momentum distribution and 
ri{k) is an uncorrelated uniform random number in the 
interval 77 G [0,27r[. In the simulations presented in Fig. 
2J the value /cq = 0.bjj.rm~^ was chosen, but we have 
checked that the results are not sensitive to this width. 
The phase of such a state in real space has a number of 
singularities that increases with the width of the state in 
momentum space. In the absence of interactions, we can 
compare the results of our numerical simulations directly 
with Q. This is shown in Figure [2] (Ic) . The numerical 
simulations are consistent with our analytic calculations. 

When interactions are turned on, we numerically ob- 
serve the expected narrowing of the momentum distribu- 
tion, corresponding to the buildup of longer range spatial 
correlations. From the real space image, it appears that 
vortices become well defined objects with a characteris- 
tic size of a few jam. This stands in stark contrast to the 
noninteracting case [compare in Fig J2] panels (Ia,b) with 



panels (IIa,b)]. As discussed above, we do not expect well 
developed vortices in the equilibrium state. We therefore 

study the interference fluctuations in the transient. In 

(2) 
our simulations we have observed that the change of g} 

in time does not change the qualitative features that we 

discuss below. 

The decay of density fluctuations in the interferrogram 
[see Fig. [2] (He)] is dramatically changed with respect to 
the noninteracting case. It now varies now much more 
with the separation a between the points that are inter- 
fered. The further the points are apart, the longer the 
correlation time of the fluctuations. To check quantita- 
tively for a ballistic scaling r ~ a, we have rescaled in 
Fig. [2] (lid) the r-axis by a. The different curves taken 
at different displacements a indeed overlap for the smaller 
distances (a < 6/im). 

Strikingly, the curves for larger separation a in Fig. 
^ lid) are above the ones for lower a. This indicates 
possible diffusive behavior of vortices at large distances 
[19] . The inset in Fig. |2]lld) shows that a diffusive scaling 
T ^ o? works indeed well for larger separation. For the 
largest distances, there is again a slight deviation toward 



larger correlation times. 

Comparing in Fig. ^the panels (I) with the panels (II), 

(2) 

it is clear that the observable gj ^ is not simply related to 
the presence of dislocations in the phase profile. Dislo- 
cations in the phase profile occur both in the interacting 
and noninteracting case, but the phase dynamics is very 
different. Our interpretation of the different phase dy- 
namics of the following. In the noninteracting case, the 
bose gas can be described in terms plane wave excita- 
tions. This is precisely the assumption that leads to the 
prediction ([2|. Vortices, i.e. zeros of the wave function, 
are incidental and bear not much physical significance. 
The situation changes when interactions are introduced. 
Plane waves are then no longer the good excitations to 
characterize the system. This is witnessed very clearly by 
the dynamics of the interference that changes dramati- 
cally between panels (I) and (II). The dynamics shown 
in panels (II) can clearly not be explained by the pre- 
diction ([2| for uncorrelated plane waves. In a real space 
picture, the correlations between the plane waves gives 
rise to the appearance of vortices with a characteristic 
size (the healingh length) as elementary excitations. 

These intensity fluctuations of the interferrogram thus 
give information that is complementary to the direct 
measurement of the condensate phase that identifies the 
zeros of the wave function. This observable could there- 
fore be of interest for ultracold atoms as well. Experimen- 
tally, it could be measured by extracting atoms from two 
spatially separated regions in the condensate and count- 
ing them with a single atom counter, as in Ref. [20 . 

It is also interesting to investigate the case of a well 
condensed quantum fluid with small fluctuations on top. 
For the numerical study, we have taken an initial con- 
dition with 96 percent of the particles in the conden- 
sate and a Gaussian distribution in momentum space for 
the thermal component. In this state with good spatial 
coherence, the fluctuations are too weak to induce free 
vortices. The Bogoliubov excitations are weakly occu- 
pied and can be considered to be a noninteracting gas. 

(2) 

We thus expect the decay of Qj ^ to depend less on a. 
Panel (IIIc) in Fig. [2] shows that the dependence on a 
is indeed much weaker than in panel (He). Especially 
when the time r is rescaled with the separation a, this 
becomes very clear. Again, as in the noninteracting case, 
the curves with a longer separation decay faster in this 
graph than the ones with a shorter separation. In this 
sense, we recover the behavior of the noninteracting gas. 
Note however that the magnitude of the interference fluc- 

(2) 

tuations is much enhanced: the g} drops by a factor of 
five, which is much larger than the factor of two that we 
obtain for the Gaussian fluctuations in the noninteract- 
ing case of panels (Ic,d). Note that for long separations, 
the curves in panel Illd) approach each other. This is 
indeed expected for phase fluctuations that have a linear 
dispersion. 

The density fluctuations of the interferrogram can be 
experimentally measured with standard quantum optical 
techniques. For this reason, exciton-polariton quantum 
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FIG. 3: The same as the corresponding panels in Fig. [2] but 
for a finite polariton life time j/{gn) — 0.01. 



fluids in semiconductor microcavities appear as the most 
natural candidate for studying the phase dynamics in this 
particular way. A specific aspect of this quantum fluid is 
the finite life time of the microcavity exciton-polaritons 
pT] , Even in the atomic case, a measurement of density 
fluctuations requires the extraction of atoms from the 
condensate, limiting the life time of the condensate as 
well. 

Thanks to progress in semiconductor manufacturing 
technology, the life time of microcavity polaritons has 
been pushed to the 100 ps regime, equivalent to a line 
width on the order of 10 /ieV, that has to be compared 
to an interaction energy that can be of the order of 1 eV. 
We thus expect that in the best microcavities, the effects 
of a finite life time can be minor. To verify this expecta- 
tion theoretically, we have repeated our simulations with 
a finite polariton life time. Numerically, this was im- 
plemented by a the inclusion of an imaginary part in the 
polariton dispersion e(k). Fig. [sjshows the results of this 
numerical simulation. The difference between the panels 
II) and III) is clear, but it is reduced with respect to the 
case of infinite life time in Fig. [2J It also appears that 
that the life time is not long enough in order to observe 
the diffusive scaling [cf. inset of Fig. |2] lid] . Actually at 
long times, the quantum fluid enters the noninteracting 
regime. This is clearly visible in panel lid), where the 
order of the curves reverses when r/a > 5ps//im. 



IV. CONCLUSIONS 

In conclusion, we have suggested to study the density 
fluctuations of the interferrogram to characterize differ- 
ent regimes of a quantum degenerate Bose gas. In partic- 
ular, a very long correlation time, that shows a ballistic 
scaling at short distances and a diffusive one at larger dis- 
tances is observed in a quantum fluid with proliferating 



vortices. These features mark a distinct contrast with 
either a noninteracting gas or a Bose-Einstein conden- 
sate. The numerical evidence for a nontrivial dynamics 
of these density fluctuations of the interferrogram urges 
for a deeper theoretical investigations. This could allow 
for a more precise characterization of different regimes in 
quantum turbulence. 
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A Fourier transform of the fields in Eq. (IT]) yields 
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where Cj = e(kj). When the phases of the different momentum components are random, the only terms that remain 
are the usual Wick contractions for which ki^2 = k4^3 and ki^2 = k3^4. The first combination gives a contribution that 
is independent of time, where the second one gives the time dependence. The time dependent part can be written as 
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Because for the non-interacting Bose gas the field '^(x) — 0) = 2g\ ^(a, t ,r -^ oo). With this relation, we obtain 
'0(x + a) has a Gaussian distribution, its density fluctu- from Eq. (A.2) the expression ^. 



(2). 



ations show bunching by a factor of two: g\ ^(a, t,r 
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